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Abstract 

Two very different methods - exact diagonaiization on finite chains and a variational method - 
are used to study the possibility of a metal-insulator transition in the symmetric half-filled periodic 
Anderson-Hubbard model. With this aim we calculate the density of doubly occupied d sites (z^) 
as a function of various parameters. In the absence of on-site Coulomb interaction (Uf) between / 
electrons, the two methods yield similar results. The double occupancy of d levels remains always 
finite just as in the one-dimensional Hubbard model. Exact diagonaiization on finite chains gives 
the same result for finite U / , while the Gutzwiller method leads to a Brinkman-Rice transition at 
a critical value (Uf), which depends on Uf and V. 
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I. INTRODUCTION 



The periodic Anderson-Hubbard model defined by the Hamiltonian 





3, <? 3 



(i) 



is meant to describe the physics of systems in which two types of electrons, one filling a 
relatively broad conduction band, the other a narrow band, are allowed to hybridize. In 
what follows we call them d and / electrons. The interactions within the bands are denoted 
by U d and Uf, respectively. In the present work we restrict ourselves to the half-filled 
paramagnetic case, that is, when there are — — N up- and down-spin electrons in an 
arbitrary dimensional lattice with N lattice sites, each of which has one d and one / orbital. 
Moreover, we restrict ourselves to the symmetric case, where an equal number of d and / 
electrons is present on the average. This is realized when Sf — (Ud — Uf)/2, if the energy is 
measured from the center of the d band [1]. 

In our previous study of this model [1] , in which we used two very different methods: a 
variational calculation using the Gutzwiller type wave function and exact diagonalization, 
we were mainly interested in the effect of the on-site interaction Ud between conduction 
electrons on the /-electron physics. In the present study we examine the effect of d-f 
hybridization (V) and of the interaction Uf between / electrons on the Hubbard physics, 
that is on the eventual metal-insulator transition at half filling. 

In the Gutzwiller-type treatment of the half-filled Hubbard model, the metal-insulator 
transition, which is known in this case as the Brinkman-Rice transition [2] , occurs at a finite 
Ud, where the number of doubly occupied d sites becomes zero. A similar transition was 
obtained by the Gutzwiller method in the half-filled periodic Anderson-Hubbard model, too 
[1], when the / electrons in the narrow band are strongly correlated. 

In contrast to that, exact diagonalization of the half-filled periodic Anderson-Hubbard 
model on finite chains gave a finite number of doubly occupied d sites for any U d , just as in 
the one-dimensional half-filled Hubbard model, where this number is finite for arbitrary U d 
[3], even though the ground state is conducting only for U d — and insulating for U d > 0. 

In this paper we will consider the half-filled symmetric Anderson-Hubbard model in 
the paramagnetic regime in the full U d > 0, Uf > sector, when it is not necessarily 
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in the strongly correlated Kondo region, and will study the possibility of metal-insulator 
transition. We should note here that in the present model a "metallic" phase is in fact a 
band insulator with hybridization gap. Therefore, we should be speaking about insulator- 
insulator transition, though its physics is the same as a metal-insulator transition due to 
the interactions between electrons. We calculate the number of doubly occupied d sites as a 
function of Ud for various values of V and Uf by both methods at the symmetric Anderson 
point, where the average number of / and d electrons per site (denoted by n/ and rid, 
respectively) is exactly 1, and examine the effects of these couplings on the one- and higher 
dimensional Hubbard physics. 

II. CALCULATION BY EXACT DIAGONALIZATION 

First, we perform exact diagonalization of the model on finite chains, where the kinetic 
energy of conduction electrons moving along the chain is described by hopping between 
nearest-neighbor d orbitals with hopping rate t. 

We consider the periodic Anderson-Hubbard model on a chain consisting of six sites with 
periodic boundary conditions. The results are shown in Figs. 1 and 2, where the density 
of doubly occupied d sites, v A , is shown as a function of U d for V = 0.1W and 0.3H^, 
respectively (W = At is the bandwidth), and for Uf — and 5W. The values calculated 
with the Bethe Ansatz for the pure Hubbard model are also shown in the figure by a solid 
line. 

One can see that when the conduction electrons of the d band are hybridized with non- 
interacting electrons in the / band (Uf = 0), the larger the hybridization the more the 
number of doubly occupied sites. The curves are reasonably close to the results obtained 
by the Gutzwiller method. The agreement gets better for stronger hybridization, while for 
weak hybridization it holds for small Ud values only. 

The values of v& decrease for finite U f and get close to those of the pure Hubbard model 
for large Uf. This suggests that the ci-electron subsystem becomes decoupled from the / 
electrons when the / electrons are strongly correlated. 

The results in this section are valid for a chain, for a one-dimensional model. In the next 
section we discuss a variational method, which might be relevant for higher dimensional 
models. 
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1 2 3 4 5 6 
U d /W 

Fig. 1. Ud as a function of Ud for V = Q.1W. The empty and filled circles indicate the results 
of exact diagonalization for Uf/W = and 5, respectively. The dashed line is the result of the 
Gutzwiller method for Uf/W = 0. The solid line is the exact solution of the one dimensional 
Hubbard-model. 
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Fig. 2. The same as Fig. 1 except that V = 0.3W. 



III. VARIATIONAL CALCULATION 

We summarize the main steps of the variational calculation following Ref . [4] . The trial 
wave function is chosen in the form 

i*> = n n + io>, (2) 

k cr 
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where the Gutzwiller projectors, Pq and Pq, which contain the variational parameters r] d 
and rjf, are written as 



^ = n t 1 " " ^"SJ > 



"r'i 



(3) 

(4) 



The variational parameters, r\ d and rjf, depend on U d and Uf, respectively. Performing the 
optimization with respect to the mixing amplitudes, and Vk, we get 



8 = n S qd£d ^ + \J[q d e d {k)~e s ] 2 + AV 2 



fceFS 



(5) 



for the ground-state energy density, where q d denotes the kinetic energy renormalization 
factor of d electrons given by 



Qd 



(1-f) 



I'd 



y - Vij v d + ^/(y - ^d) (1 - n d + ^) 



(6) 



which is formally identical to the expression found in the Hubbard model [5]. The renor- 
malized hybridization amplitude is now V = V \/qMf] the other notations are the same as 
in our previous paper [4], and the self-consistency condition is given by 



- V 

N ^ 



fceFS 



1 + 



q d e d (k) -e f 



1. 



(7) 



[q d e d (k)-e f ] + W\ 

The summation over k can be carried out assuming a constant density of states, p(e) = 1/W, 
in the interval e G [-W/2, W/2\. The values of Uf, and v d are obtained by optimizing the 
energy density with respect to these parameters numerically. For V <^W and n d = rif = 
1 the optimization condition with respect to v d and Vf results in the following coupled 
equations: 



Ud 

w 



4 \Wj q d 



U -l + 2(^ \ In 
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8(l-4i/ d ) = 0, 

(l-4i/,) = 0. 



(8) 
(9) 



When either of U d or Uf is zero, the equations are decoupled. Note that in the absence of 
hybridization Eq. (8) reduces to the result for the ordinary Hubbard model. 



We now turn to the discussion of the ci-electron subsystem. We calculate the density of 
doubly occupied d sites, u d , as a function of XJ A . The results are displayed in Fig. 3 for 
several values of Uf for a relatively weak hybridization, V = 0.1W, and in Fig. 4 for a 
stronger hybridization. 




U d /W 

Fig. 3. u d as a function of Ud- The solid, dashed, dotted and dot-dashed lines correspond to 
Uf/W = 0, 0.3, 0.5, 5 respectively. V = 0.1W in all cases. 




U d /W 

Fig. 4. v<i as a function of Ud- The solid, dashed and dotted lines correspond to Uf/W = 0, 1, 10 
respectively. V = 0.3W in all cases. 

First, we find that for Uf = 0, v d never becomes zero, just as it was obtained by the exact 
diagonalization for a finite chain, that is, the Brinkman-Rice transition does not occur. We 
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can calculate the asymptotic behavior of v d for large value of U d from the following analysis. 
When Uf — 0, qf — 1 irrespective of U d , and Eq. (8) can be solved for v d . For Ud 3> W in 
leading order we arrive at: 
2V 2 

u d oc — . (10) 

Second, we find that the Gutzwiller method leads to a Brinkman-Rice transition for any 
Uf > 0, that is, there exists a finite value, U%, where v d becomes zero. For U d > and Uf > 
the optimization conditions, Eqs. (8) and (9) for Vf and Vd are coupled, and thus both v d 
and Vf decrease when either of the interactions increases. Therefore, even for very small Uf, 
when U d is large enough, Vf also becomes small, and finally both u d and Vf simultaneously 
become zero at U d . As a matter of fact, U d and Uf play a rather similar role. If we fix U d 
at a certain value larger than 2W, a Brinkman-Rice transition occurs at a certain Uf. The 
difference in the condition (U d > 2W and Uf > 0) necessary for occurrence of a transition 
is due to the different widths of the d and / bands (W and 0, respectively) in the present 
model. 

Third, when Uf is large enough and Vf is exponentially small even for small U d , that is, 
when the system is the Kondo regime, the v d — U d curves become straight lines and coincide 
to the known behavior of the Hubbard-model. In the limit Uf 3> W we obtain: 

v d — ; r, 11 

4 8(iy + 4E K ) v ; 

where 

E K = ^J-^I^\ (12) 

K 2 l \ 16V 2 /W J V ; 

The scenario is the same for weak or strong hybridizations. 

The results described above indicate that a phase boundary can be defined in the three- 
dimensional parameter space of U d , Uf and V, which separates the region where v d — Vf — 
from that where both u d and Uf are finite. The former is a Mott insulator region and the 
latter is a hybridized band insulator region (for small Uf) or a Kondo insulator region (for 
large Uf). The phase boundaries in the U d -Uf plain are shown in Fig. 5. (About the 
boundary between a hybridized band insulator and a Kondo one, see Ref. [1].) 



IV. CONCLUSIONS 
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Fig. 5. The phase boundaries in the U$-Uf plain separating the metallic and insulating regimes. 
The solid and dashed lines correspond to V = 0.1W and 0.3W, respectively. 

In this paper we discussed the periodic Anderson-Hubbard model focusing our attention 
on the physics of the conduction electron subsystem, using two different methods: exact 
diagonalization on finite chains and a variational method of the Gutzwiller-type. We studied 
the effects of the d-f hybridization (V) and of the on-site interaction between / electrons 
(Uf) on the number of doubly occupied d sites, v&. When Uf = 0, both methods gave similar 
results. For larger Uf, however, the results of exact diagonalization in the one- dimensional 
model showed that v d approaches that obtained from the Bethe-Ansatz solution of the pure 
Hubbard model, while the Gutzwiller method indicates a Brinkman-Rice-type scenario for 
a metal-insulator transition. 

It is interesting from theoretical point of view that v& oc l/t/J according to the Bethe- 
Ansatz solution for Ud ^> W , while the Gutzwiller method gives a slower asymptotic behav- 
ior, u d oc l/Ud for Uf = 0. 
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